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^ ' Abstract 

' This paper reports the values of contributions to the electron g—2 from 300 Feynman diagrams 

. of the gauge-invariant Set Ill(a) and 450 Feynman diagrams of the gauge-invariant Set Ill(b). The 

I evaluation is carried out in two versions. Version A is to start from the sixth-order magnetic 

I anomaly Mq obtained in the previous work. The mass-independent contributions of Set Ill(a) 

^ ■ and Set ni(b) are 2.1275 (2) and 3.3271 (6) in units of (a/vr)^, respectively. Version B is based 

. on the recently-developed automatic code generation scheme. This method yields 2.1271 (3) and 

\ 3.3271 (8) in units of (a/vr)^, respectively. They are in excellent agreement with the results of 



X 



the first method within the uncertainties of numerical integration. Combining these results as 
statistically independent we obtain the best values, 2.1273 (2), and 3.3271 (5) times (a/vr)^, for 



\ the mass-independent contributions of the Set Ill(a) and Set Ill(b), respectively. We have also 

evaluated mass-dependent contributions of diagrams containing muon and/or tau-particle loop. 
Including them the total contribution of Set Ill(a) is 2.1349 (2) and that of Set Ill(b) is 3.3299 (5) 
in units of (a/vr)^. The total contributions to the muon 17—2 of various leptonic vacuum-polarization 
loops of Set Ill(a) and Set Ill(b) are 112.418 (32) and 15.407 (5) in units of (a/Tr)^, respectively 

PACS numbers: 13.40.Em,14.60.Cd,12.20.Ds,06.20.Jr 
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I. INTRODUCTION 



The anomalous magnetic moment g — 2 of the electron has played the central role in 
testing the validity of quantum electrodynamics (QED) as well as the Standard Model. The 
latest measurement of = {g — 2)/2 by the Harvard group has reached the precision of 
0.24 X 10"^ |l, 2|: 

ae(HV08) = 1 159 652 180.73 (0.28) x 10"^^ [0.24ppb] . (1) 
At__present the theoretical predic ;ion consists of QED corrections of up to the eighth order 



5|, and hadronic corrections [6 
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scaled down from 



12| and electro-weak corrections 
their contributions to the muon g—2. To compare the theory with the measurement ([1]), we 
also need the value of the fine structure constant a determined by a method independent 
of g—2 . The best value of such an a has been obtained recently from the measurement of 
h/mjQ^, the ratio of the Planck constant and the mass of Rb atom, combined with the very 



precisely known Rydberg constant and m-^\^/me'- 16| 

a"^(RblO) = 137.035 999 037 (91) [0.66ppb]. (2) 

With this a the theoretical prediction of becomes 

ae(theory) = 1 159 652 181.13 (0.11)(0.37)(0.77) x 10"^^ (3) 

where the first, second, and third uncertainties come from the calculated eighth-order QED 
term, the tenth-order estimate, and the fine structure constant ([2]), respectively. The theory 
(|3]) is thus in good agreement with the experiment ([1]): 

ae(HV08) - ae(theory) = -0.40 (0.88) x 10"^^ (4) 

proving that QED (Standard Model) is in good shape even at this very high precision. 

An alternative test of QED is to compare a(RblO) with the value of a determined from 
the experiment and theory of g—2: 

a-i(ag08) = 137.035 999 085 (12)(37)(33) [0.37ppb], (5) 

where the first, second, and third uncertainties come from the eighth-order QED term, 
the tenth-order estimate, and the measurement of ae(HV08), respectively. Although the 



uncertainty of a~^{ae08) in ([5]) is a factor 2 smaller than a(RblO), it is not a firm factor 
since it depends on the estimate of the tenth-order term, which is only a crude guess jl?! . 
For a more stringent test of QED, it is obviously necessary to calculate the actual value of 
the tenth-order term. In anticipation of this challenge we launched a systematic program 
several years ago to evaluate the complete tenth-order term 
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20|. 



The lOth-order QED contribution to the anomalous magnetic moment of an electron can 
be written as 



,(10) 



TT 



(6) 



where m^/m^ = 4.836 331 71 (12) x 10"^ and me/m^ = 2.875 64 (47) x 10"^ pT,]. In the 
rest of this article the factor (^j will be suppressed for simplicity. 

The contribution to the mass-independent term A^^^^ can be classified into six gauge- 
invariant sets, further divided into 32 gauge-invariant subsets depending on the nature of 
closed lepton loop subdiagrams. Thus far, numerical results of 27 gauge-invariant subsets. 



which consist of 3106 vertex diagrams, have been published 18|, I2ll-l26|. Five of these 27 



subsets were also calculated analytically (271, l28| . Our calculation is in good agreement with 
the analytic results. 

In this paper we report the evaluation of the tenth-order lepton g — 2 from two gauge- 
invariant subsets called Set Ill(a) and Set Ill(b). These diagrams are built from the magnetic 
moment contribution Mg (shown in Fig. [T]) which consists of 50 proper sixth-order vertices 
of three-photon-exchange type, namely diagrams without closed lepton loops (and called 
g-type. See Ref. j4| for the definition of g-type.), by insertion of various lepton vacuum- 
polarization loops. 




FIG. 1: The sixth-order q-type diagrams. The solid line represents the electron in a constant 
magnetic field. The time reversal diagrams of Mqd and Mqq are omitted for simplicity. 
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Ill(a) Ill(b) ni(c) 



FIG. 2: Typical diagrams of Set III. 

Set Ill(a). Diagrams obtained by inserting two second-order vacuum-polarization function 
n2's in Mg. The number of vertex diagrams contributing to A^^^^ is 300. 

Set Ill(b). Diagrams obtained by inserting the fourth-order vacuum-polarization function 
114 in Mg, where 114 is the sum of three fourth-order vacuum polarization loops. The 
number of vertex diagrams contributing to is 450. 

Another set ( Set III(c) of Fig. [2]) consists of diagrams obtained by inserting a light-by- 
light scattering subdiagram A4 in Mg. The total number of these diagrams contributing to 
A^^^-' is 390. Since it has a structure different from those of Sets Ill(a) and Ill(b), it will be 
treated in a separate paper. 

Evaluation of Set III (a) and Set Ill(b) is carried out in two ways. Version A is to start 
from the FORTRAN code of the sixth-order anomalous magnetic moment Mg, which was 



obtained in previous works [29|] and known to give the result identical with the analytic 
result [30|. It is thus easy to establish the validity of these FORTRAN codes for Sets Ill(a) 
and Ill(b). 

We also evaluate these sets by an alternative method, Version B, using FORTRAN 
codes generated from scratch by the recently developed automatic code generation scheme 



19| . This approach deals with the UV renormalization as well as IR subtraction terms 
as integral parts of automation. In carrying out this automation scheme, we found it useful 
to construct IR subtraction terms in a different manner from that of Version A [19]. Thus, 
Version B provides an independent confirmation of Version A. At the same time it helps 
to verify the automated code generation scheme, which is developed primarily to deal with 
the vastly more difficult problem of Set V, which consists of 6354 vertex diagrams with pure 
radiative correction. 

As is well-known, the insertion of vacuum-polarization loop such as 11*^^^ and 11*^^) in an 
internal photon line of momentum q can be expressed as a superposition of massive vector 
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ni(ad) 



ni(a,) 



FIG. 3: Typical tenth-order diagrams of Set Ill(a) obtained by insertion of two second-order 
vacuum-polarization loops 112 in lepton diagrams of the three-photon-exchange type. The subset 
Ill(a^) consists of diagrams in which 112 are inserted in different photon lines, while the subset 
Ill(as) consists of diagrams in which 112 are inserted in the same photon line. There are 150 
diagrams in each subset. 

propagators 

r (7) 

where m is the mass of the lepton forming the closed loop and a is the square of mass 
of the vector particle and p is the spectral function. This enables us to obtain Feynman- 
parametric integrals for Set Ill(b) by simply replacing the relevant photon mass squared by 
a and integrating over a. It can also be applied to diagrams of Set Ill(a) which contain 
two vacuum-polarization loops in different photon lines. This subset of Set III (a) will be 
denoted as Set Ill(ad) henceforth. 

The Set Ill(a) also contains diagrams in which two vacuum-polarization loops are inserted 
in the same photon line, which will be denoted as Set Ill(as). For these diagrams a slight 
extension of Eq. ([7]) is required. When two vacuum polarization loops are inserted in a 
photon line of momentum q, the result, omitting integrations for simplicity, is given by the 
left-hand-side of the following equation, which can be rewritten in the form on the right- 
hand-side: 

1 ^2 1 _ 1 1 /ox 

~^ •? ~^ = 2 2 • V°J 

Note that the right-hand-side is a linear combination of propagators of mass-square Oa and 
(Tb with coefficients Oaj {oa — Oh) and —Ohj {pa — o'b) ■ This enables us to write the Feynman- 
parametric integrals for the diagrams in Set Ill(as) by a simple extension of Mg integrals. 
Eq. (IH]) can be readily extended to the case in which three or more vacuum-polarization 
loops are inserted in the same photon line. 
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These adaptations require a slight modification of the numerator function V, which, for 
Mq, is given by 

5 



V = Vo + iza + z, + zJA^ , (9) 

where Zi {i = 1, ■ ■ ■ ,5), and zj {j = a, b, c) are Feynman parameters assigned to the fermion 
propagators and the photon propagators, respectively, mg and A are masses of the electron 
and photon, respectively. Ai {i = 1, ■ ■ • , 5 ) are scalar currents flowing in the fermion line 
i (see the exact definition of in Ref. 29|] ). is expressed by the Feynman parameters 
and its expression depends on the structure of a diagram. But, the expression of V in terms 
of Ai is identical for all diagrams of Mq. 

When one vacuum-polarization function is inserted in a photon line, we must replace the 
mass square of the photon in Eq. (|9]) by p(t): 

a^^pW = ^, (10) 

where m„p is the rest mass of the fermion forming the vacuum-polarization loop and the 
interval 4m^p < a < oo of Eq. ([7]) is mapped onto (0 < t < 1) for the sake of convenience. 

When two vacuum-polarization functions are inserted in the same photon line a, it follows 
from Eq. ([H]) that the denominators must be modified as follows: 

1 V, 

V V1V2 

1 ^ - Zlp^{h)p2{t2) 



i Vq - Wozlp^{h)p2{h) - Zlp^{h)p2{h){pi{h) + P2{h)) 



(11) 



2 



where 



Vi = Vo + ZaP^{ti) + (zb + ^,)A^ z = 1, 2, (12) 



for the first or second vacuum-polarization functions. 

Throughout this article we use the exact renormalized forms of 112 and 114 instead of 
intermediately renormalized forms to take advantage of the known analytic forms of their 



spectral functions [31 1. 
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II. SET III(A) 



Diagrams belonging to the Set Ill(a) are generated by inserting two second-order vacuum- 
polarization loops 112 in the photon lines of Mg. Using an identity derived from the Ward- 
Takahashi identity |l8| and time-reversal invariance and summing up all possible insertions 
of the photon spectral function reduce the number of independent integrals from 300 to 16. 
For programming purpose it is convenient to treat Set Ill(ad) and Set lll(as) separately. 



A. Set Ill(ad) 

Let MQa,P2P2 be the magnetic moment projection of the Set Ill(arf) generated from a self- 
energy diagram Mq^ =A through H) by insertion of two electron vacuum-polarization 
loops 112 in different photon lines (see Fig. [3]). The subscript P2P2 implies that two second- 
order vacuum polarization function P2's are inserted in different photon lines of the proper 
diagram Mga. To be precise MQa,P2P2 should be written as M'^^^p^],^, where the first super- 
script h refers to the open lepton line and I2 and refer to closed lepton loops. When /i, /2, 
and /s are identical so that Mqci,P2P2 is mass- independent, we omit the superscripts for sim- 
plicity. Distinction by superscript becomes necessary in Sec. Ill A 31 where mass- dependent 
terms are treated. 



as 



1. Electron g — 2: Version A 

n Version A the renormalized contribution of the diagrams of Set Ill(arf) can be written 
32| 



H 



'[Set Ill(ad)] = ^ a6a,P2P2, (13) 

a=A 

with 

aQa,P2P2 = ^MQa,P2P2 + residual renormalization terms, (14) 
where AMQa,P2P2 is the UV- and IR-finite part of MQa,P2P2 after all divergences are removed 



by intermediate renormalization by Kg and Ir operations. See Ref. 32|| for definitions of 
i^-operation and /-operation. 

When summed over all diagrams of Set Ill(ad), the UV- and IR-divergent pieces cancel 
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out and the total contribution to a^^^^ can be written as a sum of finite pieces 



H 



[Set Ill(ad) : Ver. A] = ^M^a,P2P2 



a=A 



- 3AS2,P2AM4,P2 - 3AS2AM4,P2P2 

+ A(5m4^P2(M2*,p2[/] — M2*.P2) + A5m4^P2P2 (M24/] -M2.) 

- [AS4,P2 + 2AL4,P2 - 4Afi2A52,P2]M2,P2 

- [AE4,P2P2 + 2AL4,P2P2 - 2(A52,P2)^]M2, (15) 

where the number of vertex diagrams represented by AM6a,p2P2 is 15 for a = A,B,C, 
E,F,H and 30 for a = D,G (see Fig. [1]). We use the compactified notations for the 
magnetic moment, mass renormahzation constant, wave^function renormahzation constant, 
and vertex renormahzation constants of fourth order j32|: 



AM4 


= AM4„ + AM4fe, 


A(5m4 


= ASm^a + A6m4h, 




= Ai?4a + Ai?4fe, 


AL4 


3 

= ^^(A-L4a,i + AL^ij^i 



(16) 

1=1 

where 4a and 46 refer to fourth-order diagrams with two photons crossed and uncrossed, 
respectively, and i = 1, 2, 3 refers to three consecutive lepton lines of the diagram of type 4a 
or 46. M2* is the second-order magnetic moment with a two-point vertex insertion. M2* [I] 



is the specific limit of M2* related to /-operation defined in Refs. 32|, |33|]. The subscript P2 
in Eq. (fT5|) means that a second-order vacuum-polarization function 112 is inserted in one 
of photon lines of the proper diagram in all possible ways. Similarly, P2P2 means that two 
112 's are inserted in two different photon lines in all possible ways. 

The numerical values of AM6a,p2P2 are summarized in Table [H Numerical values of 
auxiliary integrals needed to complete the renormahzation are listed in Table HTl 

Substituting the values listed in Tables HI and HTl into Eq. ( |T5l) . we obtain 

a(^°)[Set Ill(arf) : Ver. A] = 0.941 92 (7). (17) 
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TABLE I: Contributions of the Set III(a(i) diagrams to the electron g — 2 evaluated in Version 
A. Both closed loops are electron loops, ni;' is the number of Feynman diagrams represented by 
the integral. All integrals are evaluated in double precision. First 50 iterations are carried out 
using 1 X 10® sampling points per iteration. It is then evaluated using 1 x 10^ sampling points per 
iteration for 50 more iterations. 



Integral 


np 


Value (Error) 


Sampling per 


No. of 






including np 


iteration 


iterations 




15 


-0.108 564 ( 3) 


1 X 10^,1 X 10^ 


50,50 




15 


0.107 954 (11) 


1 X 10^,1 X 10^ 


50,50 


AM6C,P2P2 


15 


0.193 333 ( 6) 


1 X 10^,1 X 10^ 


50,50 


AM6D,P2P2 


30 


0.176 456 (20) 


1 X 10^,1 X 10^ 


50,50 


AM6£;,P2P2 


15 


0.142 839 (11) 


1 X 10^,1 X 10^ 


50,50 


AMqf^P2P2 


15 


0.194 882 (15) 


1 X 10^,1 X 10^ 


50,50 


AM6G,P2P2 


30 


0.542 183 (34) 


1 X 10^,1 X 10^ 


50,50 




15 


-0.190 026 (43) 


1 X 10^,1 X 10^ 


50,50 



2. Electron g—2: Version B 

The FORTRAN programs of Version B were generated by the automation code gen- 
CODEiV with slight modification. Given one-line information specifying a diagram, gen- 
CODEiV produces a set of programs for a g-type diagram of any order of the perturbation 
theory jiol. [2^. The insertion of the vacuum polarization function in a photon line is a 
trivial task requiring modification of just a few lines of the gencodeA'^ source code. The 



-ft^-operation^ method developed in Ref. 33| can be easily automated and incorporated in 
GENCODeA [19] to deal with UV divergence. IR divergence, on the other hand, is somewhat 
differently treated. 



32 



35| successfully generates the IR sub- 



The /-operation defined in the previous work 
traction terms for a g-type diagram of up to the eighth-order of the perturbation theory. 
Actually, the /-operation works even for the tenth-order case, except that the automation 
becomes tremendously complicated. This is why we sought another way to handle the IR di- 
vergence. Namely, we deviated from the strict IR power counting, on which the /-operation 
is defined, and took a more diagrammatic approach. 
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TABLE II: Auxiliary integrals for the Set Ill(arf) and Set Ill(as) with (/i/a^s) = (eee), where 
(hhh) is defined in III A[ Six lines in the middle are for Set III (a^) and bottom four lines are for 
Set III (a^). Some integrals are known exactly. Other integrals are obtained by VEGAS integration 
34|. 



Integral 


Value(Error) 


Integral 


Value (Error) 


M2 


0.5 


M2. 


1.0 


M2.[I] 


-1.0 


AB2 


0.75 


AM4 


0.030 833 612... 


A Att).^ 


1.906 340 (21) 


AL4 


0.465 024 (17) 




-0.437 094 (21) 


Mo DO 


015 687 421 


Mo. DO 

,P2 


044 077 4 (3) 




0.010 255 3 (11) 


AM4,P2 


-0.106 707 082... 


^-82, P2 


0.063 399 266... 




0.679 769 (15) 




0.200 092 (14) 


A-B4,P2 


-0.314 320 (10) 


AM4^P2P2 


-0.026 682 (2) 


ASm4^P2P2 


0.105 075 (11) 


AL4^P2P2 


0.005 481 (8) 


A-B4,P2P2 


-0.071 017 (4) 


Af2,P2:2 


0.002 558 524... 


-M^2*,P2:2 


0.008 482 (1) 


M2.,P2:2[/] 


0.032 904 (9) 


AM4,P2:2 


-0.057 587 8 (9) 


A-B2,P2:2 


0.027 902 3 (4) 


A(5m4^P2:2 


0.439 326 (81) 


AL4^P2:2 


0.094 940 (26) 


Ai?4,P2:2 


-0.199 173 (89) 



The new scheme to deal with the IR divergence, called //-R-subtraction, consists of two 
parts: One is the i?-subtraction that removes the UV-finite part of mass-renormalization 
term, which is the cause of linear IR divergence. (The UV-divergent part of the mass renor- 
malization is removed by the i^'-operation.) Once the mass renormalization is completed, 
the remaining IR divergence is only logarithmic and is easily subtracted by the second part 
called /-subtraction. This /-subtraction is similar to the previous /-operation, except that 
it uses the finite part of a vertex renormalization constant in addition to the logarithmic 
IR-divergent part as an IR-counter term. The ///2-subtraction can be readily incorporated 



in gencodeA^ 



20|. 



As far as the sixth-order diagrams are concerned, two methods of IR treatment, /- 



10 



operation or J/i?-subtraction, work fine making no significant difference. The difference 
is only finite amount in the amphtude of the magnetic moments, which can be identified an- 
alytically. Taking it into account, we obtain the relation of the magnetic moment amplitudes 
in Version A and Version B as follows: 

P2P2 ~ ,P2P2 ~ ^-^4fe,2,P2^2,P2 ~ AL4f,^2,P2P2^2 

- A(5m4fe,P2(M2.P2 - M2.P2[/]) - A5m4fe,P2P2(M2* - M24/]), 

^^6§P2P2 = ^M^C,P2P2 - ^Sm^a,P2 {M2*P2 - P2 [I] ) 

- A5m4a,P2P2(M2* -M24/]), 

P2P2 ~ P2P2 ^ 2AL4a_i P2M2,P2 — 2 AL4a^i P2P2M2 , 

P2P2 ~ ^■^6E,P2P2 A-^4a,2,P2^2,P2 AI/4a 2,P2P2^2 , 
/^IV1qpp2P2 — ^^''■'6P,P2P2' 



'-^^^^%G,P2P2 ~ '^^'^6G,P2P2' 



^^m,P2P2 - ^^QH,P2P2- (18) 

Note that the Version B of AMQa,P2P2 absorbs not only A6m terms but also part of AL4 
terms. From Eqs. (fT5|) and (fT8|) we obtain 



H 

^(i°)[Set lll(a,) : Ver. B] = AM^^^ 

a=A 

- 3A52,P2AM4,P2 - 3A52AM4,P2P2 

- [A54,P2 + AL4,P2 - 4A52A52,P2]M2,P2 

- (A54,P2P2 + AL4,P2P2 - 2(A52,P2)^)M2. (19) 

Of course this shift of terms in Eq. (fTSll does not affect the final result. 

This is a trivial change for the Set 111. However, in Set V, which consists entirely of 
g-type tenth-order diagrams, residual renormalization terms of A6m type give rise to linear 
IR-divergences which complicate the analysis of the renormalization scheme. Thus there is 
an advantage in removing the self-mass terms completely, not just their UV-divergent parts. 

Substituting the values listed in Tables HTl and lllll in Eq. (fT9l) . we obtain 

a(^°)[Set lll(arf) : Ver. B] = 0.941 81 (12), (20) 
which is in good agreement with ( 1T7|1 . 
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TABLE III: Version B contributions of the Set Ill(ad) diagrams to the electron g — 2. The 
corresponding programs are created by GENCODEiV. Both closed loops are electron loops, np is the 
number of Feynman diagrams represented by the integral. All integrals are evaluated in double 
precision. 



Integral 


np 


Value (Error) 


Sampling per 


No. of 






including np 


iteration 


iterations 




15 


-0.130 613 (29) 


1 X 10^ 


700 


^MqB,P2P2 


15 


-0.031 263 (26) 


1 X 10^ 


700 


^Mqc,P2P2 


15 


0.080 981 (24) 


1 X 10^ 


700 




30 


0.170 496 (41) 


1 X 10^ 


700 


^Mqe,P2P2 


15 


0.183 485 (29) 


1 X 10^ 


700 


AMeF,P2P2 


15 


0.194 756 (31) 


1 X 10^ 


700 


^Mqg,P2P2 


30 


0.541 900 (69) 


1 X 10^ 


700 


^Mqh,P2P2 


15 


-0.189 816 (54) 


1 X 10^ 


700 



3. Mass-dependent terms A2 and of Set III (ad) 

Once FORTRAN programs for mass-independent Set Ill(ad) diagrams are obtained, 
it is straightforward to evaluate contributions of mass-dependent term A2^°''(me/m^), 
A''2^\me/mr), and A^r^'^\'me/m^,me/mr). We just have to choose an appropriate fermion 
mass rrLyp in Eq. ( [T0|) . Obviously the residual renormalization terms of Set III(a(i) are slightly 
more complicated because of insertions of vacuum-polarization loops in several photon lines. 
Note also that the integrands of these sets may be strongly peaked because of their depen- 
dence on {me/m^y or [me/mT-Y which makes them more susceptible to the digit deficiency 
problem. 

Of course we can evaluate them by either Version A or Version B. Since we have estab- 
lished their equivalence, we may choose either one, say Version A. 

In the general case [hhh), where I2 7^ ^3, the residual renormalization terms of Set Ill(arf) 
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in Version A have the form 

H 



a 



6a,P2P2 
a=A 

+ A5mlfi{Ml%[I] - mI%) + A5mlf^{Ml%[I] - mI%) 
+ A5ml%'t^^{M2^[I]-M2^) 

- [AbI^^:'^ + 2ALf;|) - AAB2AB^if^]M^^^^'^ 

- [ABl^l^i + 2ALlfi - AAB2ABlf^]M!^^i'^ 

- [ABff^^l + 2ALlf^j}, - A{AB^i^;.^^ AB^^t)]M,. (21) 

For instance, for (lihh) = (eem), the first, second, and third symbols refer to the open 
electron hne, electron loop, and muon loop, respectively. Some superscripts are denoted as 
{hh) or (/1/3) since they have only one internal loop. Superscripts (/i) on AB2, M2, etc., are 
omitted for simplicity since these terms are mass- independent. Note also that the second 
and third loops appear interchangeably in the case of Set Ill(a). Thus Eq. fl2T]) represents 
the sum of (eem) and {erne). 

If I2 and Z3 represent identical particles, duplicate terms of Eq. (l2Til must be dropped to 
avoid double counting. 

Substituting the values listed in Tables HVl and IXll into Eq. fl2T|) . we obtain 

41°) [Set III(ad)(^^"^)] = 0.003 166 (4). (22) 

The contributions of diagrams of (eet), (emm), etc., can be calculate by just changing 
the mass parameters in FORTRAN programs. The residual renormalization can be carried 
out using Eq. fl2Tl) paying attention to whether I2 = I3 or not. We present the final results 
without giving details: 

4^°) [Set III(ad)('"*)] = 0.000 031 55 (8), (23) 
a(^°)[Set Ill(ad) = 0.000 080 45 (11), (24) 
4^°) [Set III(arf)("™*)] = 0.000 005 56 (2). (25) 

These results have been confirmed by comparison with the results of Version B. 

The trend of mass dependence of these results indicates clearly that (ett) case will be an 
order of magnitude smaller than (l25l) . Thus it may be ignored at present. 
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TABLE IV: Mass-dependent contributions of the Set Ill(arf) diagrams to the electron g—2 evaluated 
in Version A. One closed loop is electron loop and the other is muon loop. Each integral is the sum 
of (eem) and (erne) types, np is the number of Feynman diagrams represented by the integral. 
All integrals are evaluated in double precision. 



Integral np Value (Error) Sampling per No. of 

including np iteration iterations 

^^4a!pIp2 30 -0.000 324 ( 1) 1 x lO'^ 20 

AMgp^P2 30 0.000 357 ( 2) 1 x lO'^ 20 

AMgpJp2 30 0.000 496 ( 1) 1 x lO'^ 20 

AM^'j^'^lpr, 60 0.000 533 ( 3) 1 x lO'^ 20 

AMgp^P2 30 0.000 316 ( 2) 1 x lO'^ 20 

AM^pJilp^ 30 0.000 659 ( 3) 1 x lO'^ 40 

^^6G,P2P2 60 0.001 657 (12) 1 x lO'^ 60 

AM(5^p^2P2 30 -0.000 313 (12) 1 x lO'^ 20 



4. Muon g-2. Set Ill(ad) 

The leading contribution to the muon g — 2 comes from the (mee) case where both loops 
consist of electrons, and m stands for the open muon line. Results of numerical evaluation 
in Version A are listed in Table |Vl From this Table and Table I VI I we obtain 



a 



.{,^°)[Set III(arf)("^"^)] = 42.460 4 (188). (26) 



Next largest contribution comes from [mme). We list only the result: 



,(10) f 



a-"^[Set III(ad)(™™^)] = 11.416 9 (22). (27) 



We also obtained 



a},i"^[Set Ill(ad)^"'*^] = 0.421 97 (19), (28) 
a(,^°)[Set III(ad)('"™*)] = 0.110 71 (1), (29) 
4^°) [Set Ill(ad) = 0.008 50 (10). (30) 



These results are in good agreement with the results of Version B. 
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TABLE V: Contributions of the Set Ill(arf) diagrams to the muon g — 2 evaluated in Version A. 
Both closed loops are electron loops, ni;' is the number of Feynman diagrams represented by the 
integral. All integrals are evaluated using real(lO) arithmetic built in gfortran. 



Integral 




Value (Error) 


Sampling per 


No. of 






including np 


iteration 


iterations 



, ^(mee) 


15 


-35.558 8 (48) 


1 X 10^,4 X 10^ 


20, 


180 


^^'^6S,P2P2 


15 


44.427 6 (68) 


1 X 10^,4 X 10^ 


20, 


180 


. (mee) 
^^'^6C,P2P2 


15 


19.208 7 (65) 


1 X 10^,4 X 10^ 


20, 


180 


^^'^6D,P2P2 


30 


28.117 1 (70) 


1 X 10^,4 X 10^ 


20, 


280 


^^'^6£;,P2P2 


15 


30.980 7 (64) 


1 X 10^,4 X 10^ 


20, 


80 




15 


18.790 9 (60) 


1 X 10'^,4 X 10^ 


20, 


180 


^^'^6g,P2P2 


30 


58.890 (71) 


1 X 10'^,4 X 10^ 


20, 


180 


^^'^6/i,P2P2 


15 


-51.550 (70) 


1 X 10^,4 X 10^ 


20, 


280 



TABLE VI: Auxiliary integrals for Set Ill(arf) and Set Ill(as), where {hhh) - 
integrals are known exactly. Other integrals are obtained by VEGAS integration 



(mee). Some 



Integral 



Value(Error) 



Integral 



Mtpl 


1.094 275 (44) 


^(me) 
^'^2*,P2 




-2.183 21 (16) 


. „(me) 
^-"2,P2 




11.151 07 (49) 


A 

LJ^lVl^p2 


A r 

^-'^4,P2 


3.119 86 (66) 




. (mee) 
LAiW^ P2P2 


-1.959 37 (30) 




A r (™ee) 
^-^4,P2P2 


4.960 40 (63) 


A nimee) 
^-°4,P2P2 



Value(Error) 



2.349 75 (29) 
1.885 733 (16) 
-0.628 831 80 (2) 
-3.427 88 (49) 
16.575 79 (52) 
-6.353 75 (62) 



M 
M 
AM, 
AL 



(mee) 
2,P2:2 
(mee) r^-, 
2*,P2:2l-'J 
(mee) 
4,P2:2 
(mee) 
4,P2;2 



2.718 651 (90) 
-5.107 35 (28) 
-3.484 52 (83) 
10.621 5 (13) 



(mee) 
2*,P2:2 
(mee) 
2,P2:2 
(mee) 
4,P2:2 
(mee) 
4,P2:2 



M. 
AB. 
A5m 
AB 



6 


162 


33 (39) 


5 


330 


35 (12) 


35 


742 


2 (12) 


-12 


811 


9 (12) 
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B. Set III(a,) 



1. Electron g — 2. Version A 

Let MQa,P2:2 be the magnetic moment projection of the Set Ill(as) generated from self- 
energy-hke diagrams 6a {a =A through H) by insertion of two 112 's in the same photon hne 
(see Fig. |3]) . The renormahzed contribution due to these diagrams can be written in a way 
similar to Eq. f|T3|) . 

When summed over all the diagrams of Set Ill(as), the UV- and IR-divergent pieces 
cancel out and the total contribution to a^^°^ can be written in Version ^4 as a sum of finite 



pieces (which is similar to Eq. (5.39) of Ref. |32|): 



H 

a(^°)[Set III(a,) : Ver. A] = J] AM6a,P2:2 

a=A 

- 3A52,P2:2AM4 - 3Afi2AM4,P2:2 

+ A(5m4(M2.,P2:2[/] - M2*,P2:2) + A5m4,P2:2 (M24/] - M2* ) 

- [A54 + 2AL4 - 2(A52)']M2,P2:2 

- [A54,P2;2 + 2AL4,P2;2 - 4Afi2Afi2,P2:2]M2. (31) 

The numerical values of AM6q,,p2:2 are summarized in Table I VI II Numerical values of aux- 
iliary integrals needed to complete the renormalization are listed in Table [Tll 
Substituting the values listed in Tables IITI and Ivnl into Eq. (!3T!) . we obtain 

4^°) [Set III(a,) : Ver. A] = 1.185 56 (20). (32) 



2. Electron g—2. Version B 

For the reason discussed in Sec. Ill A 21 we obtain in Version B a formula for 
ai^''''[Set Ill(as)] which is different from (13T!) : 

H 

41°) [Set III(a,) : Ver. B] = AM^f ^^2:2 

a=A 

- 3A52,P2:2AM4 - 3A52AM4,P2:2 

- [A£4 + AL4 - 2(A52)=']M2,P2:2 

- (A54,P2:2 + AL4,P2:2 " 4A52,P2:2A52)M2, (33) 
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TABLE VII: Version A contributions of Set Ill(as) diagrams to the electron g — 2. np is the 
number of Feynman diagrams represented by the integral. All integrals are evaluated in double 
precision. 



Integral 


np 


Value (Error) 


Sampling per 


No. of 






including np 


iteration 


iterations 




15 


-0.204 682 (54) 


1 X 10^,1 X 10^ 


50,50 




15 


0.413 110 (56) 


1 X 10^,1 X 10^ 


50,50 


AM6C,P2:2 


15 


0.458 938 (53) 


1 X 10^,1 X 10^ 


50,50 


^Mqd,P2:2 


30 


0.281 276 (54) 


1 X 10^,1 X 10^ 


50,50 


^Mqe,P2:2 


15 


0.220 637 (23) 


1 X 10^,1 X 10^ 


50,50 


AMqf,P2:2 


15 


0.317 657 (48) 


1 X 10^,1 X 10^ 


50,50 


AMqg,P2:2 


30 


0.765 073 (87) 


1 X 10^,1 X 10^ 


50,50 


^Mqh,P2:2 


15 


-0.409 439 (98) 


1 X 10^,1 X 10^ 


50,50 



where 



^■^6A P2:2 ~ -^^eA P2:2 ~ 2^-^46,1 ^2,P2:2 ^ 2 AL4;, i^P2: 2M2, 

^^6b!p2:2 — '^^6p!p2:2 ~ ^-^4b,2^2,P2:2 " AL46_2,P2:2^2 

- A6mib{M2*P2:2 - M2*P2:2[I]) " A5m46,P2:2 (M2. - M2*[/]), 
^^6C,P2:2 = AM^gp2^2- A'^"^4a(M2.p2:2-M2.p2:2[/]) 

- A5m4a,P2;2(M2. -M24/]), 

^■^6D P2:2 ~ ^^QD P2:2 ~ 2 AL4a iM2,P2:2 2AL4a,l,P2: 2M2, 

^'^6£;!p2:2 ~ ^^6£;!p2:2 ~ ^-^4a,2^2,P2:2 ~ ^-^4a,2,P2:2^2; 

AAiWgpp2:2 ~ ^^^-'6P,P2:2' 

^^'^6G,P2:2 ~ ^^'^6G,P2:2' 

AM(gp2:2 = AM(Jp2:2- (34) 

From Tables [III and IVIIII we obtain 

4^°) [Set III (a,) : Ver. 5] = 1.185 26 (24), (35) 
in good agreement with fl32|) . 
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TABLE VIII: Version B contributions of Set Ill(as) diagrams to the electron g—2. Programs are 
created by gencodeA^. nj? is the number of Feynman diagrams represented by the integral. All 
integrals are evaluated in double precision. 



Integral 


np 


Value (Error) 


Sampling per 


No. of 






including np 


iteration 


iterations 




15 


-0.280 605 (31) 


1 X 10^ 


600 




15 


-0.068 830 (70) 


1 X 10^ 


600 


AM6C,P2:2 


15 


0.070 002 (67) 


1 X 10^ 


600 


AM6D,P2:2 


30 


0.269 937 (106) 


1 X 10^ 


600 




15 


0.297 943 (60) 


1 X 10^ 


600 


AMqf,P2:2 


15 


0.317 432 (69) 


1 X 10® 


600 


AMqg,P2:2 


30 


0.764 711 (127) 


1 X 10® 


600 


^Mqh,P2:2 


15 


-0.409 110 (103) 


1 X 10® 


600 



3. Mass- dependent terms A2 and A3 of Set Ill(as) 
For the Set Ill(as) we have (in Version A) 



41°) [Set III(a.)('^'^'«)] = Am£;S1, 

a=A 



+ ASm,{M2^,P2:2[lf''''''^ - Ml^'fil) + A<5mi'S(M2.[/] - M2.) 

- [Afi4 + 2AL4 - 2{AB2f]M^i^pij 

- [A^i^^g^ + 2/^1^^ - AAB2ABi'fi^%)M2 . (36) 
Substituting the values hsted in Tables HXl and IXll into Eq. ( l36l) . we obtain 

41°) [Set III(a,)(^'=™)] = 0.004 12 (10). (37) 



We also obtained 



a(^°^[Set III(a,)("™"^)] = 0.000 144 7 (9), (38) 
41°) [Set III(a,)(""*)] = 0.000 045 95 (23), (39) 
4^°) [Set III (a, )(""*)] = 0.000 009 88 (8). (40) 
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TABLE IX: Contributions of diagrams of Set Ill(as) containing one electron vacuum-polarization 
loop and one muon vacuum-polarization loop, np is the number of Feynman diagrams represented 
by the integral. All integrals are evaluated in double precision. 



Integral 


np 


Value (iiirrorj 
including np 


bampiing per 
iteration 


JNo. or 
iterations 


'-^^^^C,C.P2:2 
^-^^■^6E,P2:2 

^-'"6i/,P2:2 


30 
30 
30 

60 
30 
30 
60 
30 


-0.000 541 ( 4) 
0.001 434 (17) 
0.001 334 (26) 
0.000 747 (26) 
0.000 494 (10) 
0.001 233 (16) 
0.002 397 (78) 

-0.000 995 (43) 


1 X 10^ 
1 X 10^ 
1 X 10^ 

1 X 10'^ 
1 X 10^ 
1 X lO'^ 
1 X lO'^ 
1 X 10^ 


20 
20 
20 
20 
20 
20 
20 
20 



TABLE X: Contributions of the Set Ill(as) diagrams to the muon g — 2 evaluated in Version A. 
Both closed loops are electron loops, ni;' is the number of Feynman diagrams represented by the 
integral. All integrals are evaluated using real(lO) arithmetic built in gfortran. 



Integral np Value (Error) Sampling per No. of 

including np iteration iterations 



^^^^&A,P2P2 


15 


-38.157 8 (106) 


1 X 10^4 X 


10^ 


20, 80 




15 


51.107 5 (103) 


1 X 10^4 X 


10^ 


20, 180 


A ly^irnee) 
^^'^6C,P2P2 


15 


20.547 7 (113) 


1 X 10^,4 X 


10^ 


20, 180 


^^'^6D,P2P2 


30 


31.469 8 (113) 


1 X 10^,4 X 


10^ 


20, 280 


^^^^6E,P2P2 


15 


32.799 (89) 


1 X 10^4 X 


10^ 


20, 80 


. „(mee) 


15 


19.014 4 (96) 


1 X 10^4 X 10^ 


20, 180 


^-'"6G,P2P2 


30 


61.519 4 (101) 


1 X 10'^, 4 X 


10^ 


20, 180 


^-"^6H,P2P2 


15 


-55.142 4 (95) 


1 X 10^,4 X 


10^ 


20, 280 



These results are in good agreement with those of Version B. The contribution of the {ett) 
term is negligibly small. 
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TABLE XI: Auxiliary integrals which depend on the mass ratio m^/me- Those for Set III(a(i) are 
listed on the left side. Those for Set Ill(as) are listed on the top half of the right side. Those for 
Set ni(b) are listed on the bottom half of the right side. 



llllc^l cLl 




1 n 1 oo'y" o 1 
liiLtJgi cll 








— U.UUU Ulo iy y 1 J 


^^^^4,P2:2 


n nnn n7A 7 


y^i 




—U.UUU uiy i ( i) 


iW2*,P2:2 


n nnn nnn o 
U.UUU uuu y 


(r\\ 




u.uio oyu u ID j 


1V12* ^P2:2['l r 


n nnn n^Q 7 
U.UUU uoy 1 






n 0/1/1 nsQ /I n ni ^ 

U.U41:4 Uoy 41: I^IUI j 


A nCeem) 
^-"2,P2:2 


n nnn n'^fi 9 

U.UUU UoD z 




A/To no r ri C*^™-) 


u.uiu z / 41: z yzoDj 


A R(eem) 
^-"4,P2:2 


n nnn Ano 7 

— U.UUU ^\Ja i 






0.000 009 4 (00) 


A r (een^) 
^-^4,P2:2 


0.000 016 


( 59) 




-0.000 091 5 ( 4) 


^0^4 P2.2 


0.001 261 2 


(63) 


A r (e'") 


0.000 012 7 ( 6) 








A D(eem) 


-0.000 129 ( 3) 


-0.000 068 2 


(7) 


A r (eem) 
^-^4,P2P2 


-0.000 113 6 ( 6) 


mJ?4 


0.000 005 9 


(0) 




0.000 253 9 ( 5) 


Ma* ,P4 


0.000 052 


(1) 


. n- (eem) 


0.000 195 1 ( 3) 


Ai^4P4 


-0.000 322 


(9) 






A T (em) 
^-^4,P4 


0.000 054 2 


(13) 






AJm^^p] 


0.000 878 6 


(15) 



4. Muon g-2. Set Ill(sis) 

The leading contribution to the muon g—2 comes from the case where both loops consist 
of electrons, namely the (mee) case, where m stands for the muon. Results of numerical 
evaluation in Version A are listed in Table |Xl From this Table and Table I VII we obtain 

a(^°)[Set Il^a,)^"^"^)] = 43.048 8 (194). (41) 



Next leading term is 



We also have 



a(,^°)[Set III(a,)('"'^'")] = 12.190 3 (176). (42) 

a(,i°)[Set III(a,)(™^*)] = 0.469 43 (39), (43) 

a(,^°)[Set III(a,)(™"*)] = 0.150 11 (21), (44) 

a{,^°)[Set III(a,)(™")] = 0.013 68 (3). (45) 
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FIG. 4: Typical tenth-order diagrams of Set ni(b) obtained by insertion of a fourth-order vacuum- 
polarization loop 114 in lepton diagrams of the three-photon-exchange type. Altogether there are 
450 diagrams of this type. 

These results are in good agreement with those of Version B. 



III. SET III(B) 

Diagrams belonging to this set are generated by inserting a proper fourth-order vacuum- 
polarization loop 114 (consisting of three diagrams) in the photon lines of Mg. Time-reversal 
invariance and use of the photon spectral function reduce the number of independent 
integrals from 450 to 8. These integrals are represented by the "self-energy-like" diagrams 
of Fig. [H A typical diagram is shown in Fig. HI 



1. Electron g—2: Version A 

Let MQa,P4 be the magnetic moment projection of the set of diagrams generated from a 
self-energy diagram a (=A through H) of Fig. [T] by insertion of 114 and an external vertex. 
The renormalized contribution due to the Set Ill(b) diagrams can then be written as 

H 

a(^°)[Set Ill(b) : Ver. A] = (^ea,P4, (46) 

a=A 

with 

a6o,P4 = AM6q,,p4 + residual renormalization terms, (47) 
where all divergences, except those within 114, are removed by intermediate renormalization 



by Ks and operations. (See Ref. |32|.) 

The numerical values of Set Ill(b) integrals are summarized in Table Kill Numerical 
values of auxiliary integrals needed to complete the renormalization are listed in Table IXIIII 
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TABLE XII: Version A contributions of Set Ill(b) diagrams to the electron g—2. np is the number 
of Feynman diagrams represented by the integral. All integrals are evaluated in double precision. 



Integral 


np 


Value (Error) 


Sampling per 


No. of 






including np 


iteration 


iterations 


AM6A,P4 


15 


-1.275 23 ( 8) 


1 X 10^, 1 X 10^ 


50,50 


AMeB,P4 


15 


1.865 05 (14) 


1 X 10^, 1 X 10^ 


50,50 




15 


1.593 72 (14) 


1 X 10^, 1 X 10^ 


50,50 


AMeD,P4 


30 


1.166 99 (14) 


1 X 10^, 1 X 10^ 


50,50 


AMge,P4 


15 


1.212 50 ( 6) 


1 X 10^, 1 X 10^ 


50,50 


AMsF,PA 


15 


1.113 25 (14) 


1 X 10^, 1 X 10^ 


50,50 


AMgg,P4 


30 


2.948 70 (24) 


1 X 10^, 1 X 10^ 


50,50 


AMgh,P4 


15 


-2.231 76 (22) 


1 X 10^, 1 X 10^ 


50,50 



When summed over all the diagrams of Set Ill(b), the UV- and IR-divergent pieces cancel 
out and the total contribution to a^^°^ can be written as a sum of finite pieces (which is similar 



to Eq. (5.39) of Ref. [3^): 



H 

a(^°)[Set Ill(b) : Ver. A] = ^ AM6„,p4 

a=A 

- 3A52,P4AM4 - 3A52AM4,P4 

+ A5m4(M2*,P4[/] - M2.,P4) + A5m4,p4(M24/] - M2O 

- [A54 + 2AL4 - 2(Afi2)lM2,P4 

- [A54,P4 + 2AL4,P4 - 4A52A52,P4]M2. (48) 

Terms with suffix P4 in Eq. ( HHj) are obtained by insertion of 114 in the photon lines of 
diagrams. Note that i^-operation is not applied to 114 so that we have M2,p4, instead of 
AM2,P4, in Eq. 

Substituting the values listed in Tables IXiT] and IXIIII into Eq. fj48|) . we obtain 

a(^°)[Set Ill(b) : Ver. A] = 3.327 14 (56). (49) 
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TABLE XIII: Auxiliary integrals for the Set Ill(b). Some integrals are known exactly. Other 
integrals are obtained by VEGAS integration. 



Integral Value(Error) Integral Value(Error) 

M2 0.5 M2,P4 0.052 870 652... 

M2* 1.0 M2*,P4 0.145 597 (21) 

M2*[I] -1.0 M2*,P4[I] -0.016 526 (69) 

AM4 0.030 833 612... AM4,P4 -0.288 997 (12) 

ASrriA 1.906 340 (21) A5m4,P4 1-773 79 (26) 

AB2 0.75 AS2,P4 0.183 666 8 (18) 

AB4 -0.437 094 (21) AS4,P4 -0.816 23 (25) 

AL4 0.465 024 (17) AL4,P4 0.559 72 (25) 



2. Electron g—2: Version B 

Let us now treat Set Ill(b) by the method based on the automated code generation 
scheme. In this approach, the contribution from Set ni(b) is expressed as 

H 

a(^°)[Set Ill(b) : Ver. 5] = ^ AMg^,, 

cy.=A 

- 3AS2AM4,P4 - 3AS2,P4AM4 

- (A54 + AL4 - 2(A52)^)M2,P4 

- (A54,P4 + AL4,P4 - 4A52A52,P4)M2, (50) 
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where 

^^6AP4 = ^M^A,P4 - 2AL4b,lM2,P4 - 2AL4fe,l,P4M2, 

- A5m4b(M2.p4 - M2*P4[/]) - A(5m46,P4(M2. - M^^I]), 

^M^C,P4 = AM^^_^p4- A5m4a(M2.p4-M2.p4[/]) 

- A(5m4a,P4(M2* -M2.[/]), 

^^6d1p4 = ^^6d1p4 - 2AL4a,lM2,P4 - 2AL4a,l,P4M2, 
^^62p4 = AMg^gp4 - AL46,2M2,P4 - AL4fe,2,P4M2, 

/\Mg]p, = AM^^Jp,, 

AMgp, = AM^Jp,. (51) 

Using the code generator we obtained the programs of the magnetic moments M6q,^p4, 
a = A, . . . , H , and M^^, M4a,p4,, a = A,B. The programs for the renormahzation constants 
L4Q P4, L^a, B^a Pi, B^a, 5"^4o,P45 Sm^a are also automatically generated. Other quantities, 
A-B2,P4, ^2,P4 are very simple so that they are calculated by using hand- written programs. 
The values of Ai?2 and M2 are analytically known. 

The results of numerical integration by VEGAS are shown in Table IXIVi 
Substituting the numbers shown in Tables IXIVI and I XIIII into Eq. ( l50l) . we obtain 

4^°) [Set Ill(b) : Ver. B] = 3.327 07 (78) (52) 

in good agreement with (H9|) . where the uncertainty is from the numerical integration only. 
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TABLE XIV: Version B contributions of Set ni(b) diagrams to the electron g — 2. Programs are 
created by gencodeA^. np \s the number of Feynman diagrams represented by the integral. All 



integrals are evaluated 


in double 


precision. 






Integral 


np 


Value (Error) 


Sampling per 


No. of 






including np 


iteration 


iterations 


AM6A,P4 


15 


-1.673 94 (34) 


1 X 10^, 1 X 10^ 


50, 50 


AM6B,P4 


15 


-0.959 45 (24) 


1 X 10^, 1 X 10^ 


50, 50 


AM6C,P4 


15 


0.427 34 (21) 


1 X 10^, 1 X 10^ 


50, 50 


AM6D,P4 


30 


1.110 10 (32) 


1 X 10^, 1 X 10^ 


50, 50 


AM6B,P4 


15 


1.498 03 (18) 


1 X 10^, 1 X 10^ 


50, 50 


AM6P,P4 


15 


1.113 12 (20) 


1 X 10^, 1 X 10^ 


50, 50 


AM6G,P4 


30 


2.947 48 (36) 


1 X 10^, 1 X 10^ 


50, 50 


AM6H,P4 


15 


-2.231 66 (28) 


1 X 10^, 1 X 10^ 


50, 50 



3. Mass- dependent terms A2 of Set Illfb) 

The residual renormalization scheme (in Version A) for the (em) term is the following: 
4^°) [Set III(b)(^")] = f;AM(:-;, 

a=A 

+ A5m4(M2*,P4[/]("'") - M^^t%) + ASm[yl{M2^[I] - M2O 

- [A54 + 2AL4 - 2{AB2f]M^'^pl 

- [Ai?S + 2AL[yl - 4Ai?2Ai?(7i]M2. (53) 
Substituting the values hsted in Tables IXll and IXVl into Eq. (!53|) . we obtain 

41°) [Set III(b)('^™)] = 0.002 794 (1). (54) 

It is easy to obtain the contribution of tau lepton loop instead of the muon loop. We 
have simply to replace the muon mass by the tau mass in the FORTRAN programs. For 
instance a crude calculation in Version A yields 

4^°) [Set III(b)("*)] = 0.000 021 42 (1), (55) 

which is two orders of magnitude smaller than (E 
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TABLE XV: Contributions of (em) diagrams of Fig. [2] to the Set Ill(b). np is the number of 
Feynman diagrams represented by the integral. Ah integrals are evaluated in double precision. 



Integral np Value (Error) Sampling per No. of 

including np iteration iterations 

AmJ554 15 -0.000 389 ( 2) 1 x 10^ 20 

^^6b!p4 15 0.000 943 ( 6) 1 x lO'^ 20 

^^6c!p4 15 0.000 925 ( 6) 1 x lO'^ 20 

^^6d!p4 ^0 O-OOO 592 (4) 1 X lO'^ 20 

^^GE^L 15 0.000 371 ( 4) 1 X lO'^ 20 

AM^'jTp^ 15 0.000 773 ( 7) 1 x lO'^ 40 

^^GG^k 30 0.001 626 (26) 1 x lO'^ 60 

^^ITpa 15 -0.000 751 (20) 1 x lO'^ 20 



4. Muon g-2. Set Ill(b) 

The leading contribution to the muon g — 2 comes from the case containing an electron 
loop, namely the (me) case, where m stands for the muon. Results of numerical evaluation 
( Version A) are listed in Table IXVIi From this Table and Table IXVIII we obtain 

a^i°)[Set III(b)('"")] = 11.936 7 (45). (56) 

We also obtained ( Version A) 

a^i°)[Set III(b)('"*)] = 0.143 60 (1). (57) 

These results are confirmed by Version B calculation. 



IV. DISCUSSION 



As was noted earlier Version A and Version B differ in the treatment of self-energy 
subtraction and IR divergence. Furthermore, the actual algebraic form of integrands in the 
first method |29| is quite different from the second one because "Kirchhoff 's laws" satisfied by 



the scalar currents 



32] were used extensively to make the integrand as compact as possible 



to save the computing time. Thus the two calculations can be regarded as independent 



26 



TABLE XVI: Contributions of (me)-type diagrams of Set Ill(b) to the muon g — 2. np is the 
number of Feynman diagrams represented by the integral. All integrals are evaluated using the 
real(lO) arithmetic built in gfortran. 



Integral 


np 


Value (iLirrorj 
including np 


bamphng per 
iteration 


JNo. or 
iterations 


A 

A^(™e) 
^^'^6E,P4 
A 

A /i^Cme) 
^^'^6G,P4 
A 

^-'^-'6H,P4 


15 
15 
15 
30 
15 
15 
30 
15 


-16.0005 (53) 
25.3670 (61) 
2.8871 (61) 
12.8166 (79) 
13.6292 (44) 
6.9195 (59) 
24.6728 (59) 
-23.0066 (57) 


1 X 10^ 

1 X 10^ 
1 X 10^ 

1 X 10^ 
1 X 10"^ 
1 X 10^ 
1 X 10^ 
1 X 10^ 


600 

900 

900 

1000 

600 

600 

700 

800 



TABLE XVIL Auxiliary integrals for the Set Ill(b)^"*^^. Some integrals are known exactly. Other 
integrals are obtained by VEGAS integration. 



Integral 


Value (Error) 


Integral 


Value(Error) 


M2 


0.5 


p4 


1.493 651 (84) 


M2* 


1.0 




3.122 88 (16) 


M2.[I] 


-1.0 




-2.996 76 (24) 


AM4 


0.030 833 612... 




-0.438 76 (26) 




1.906 340 (21) 


A6m'(^pl 


13.651 22 (80) 


AB2 


0.75 


. „(me) 
^-°2,P4 


2.439 109 (53) 


ABi 


-0.437 094 (21) 


^-°4,P4 


-3.826 32 (71) 


AL4 


0.465 024 (17) 


^-^4,P4 


3.653 31 (42) 



of each other and the results agree within their error bars. Thus they may be combined 
statistically to yield the values listed below: 

4^°) [Set Ill(ad)] = 0.9419 (1), (58) 

4^°)[Set III(a,)] = 1.1854 (2), (59) 

a(i°)[Set 111(b)] = 3.3271 (5). (60) 
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The mass-dependent contribution of Set lll^ad) to the electron g — 2, the sum of (l22l) . 
(I23D, ([24]), and (ESD, is given by 

4^°) [Set III(ad)(mass-dep)] = 0.003 28 (1), (61) 

while the mass-dependent contributions of Set Il^a^) to is the sum of (jSZl), ([2SD, 
and (HDD: 

4^°) [Set III(a,)(mass-dep)] = 0.004 32 (10). (62) 
The total contribution of Set Ill(a) to Oe is the sum of (!58|) . (!59l) . fl6T]) . (l62l) : 

a(i°)[Set Ill(a) (all terms)] = 2.1349 (2). (63) 

Similarly, from (|5^ . fl55|) . and ( 160|) we obtain 

4^°) [Set Ill(b) (all terms)] = 3.3299 (5). (64) 

The total contribution of Set Ill(ad) to the muon g-2, the sum of ([26]), ([27]), ([28]), ([29]), 
([50]). and ([SS]), is 

aJ,^°)[Set III(ad)(all terms)] = 55.360 (19), (65) 

while the total contribution of Set III(aJ to the muon g-2, the sum of ([H]), ([52]), ([13]), (HI) . 
([15]). and ([59]), is 

aJ,^°)[Set III(a,)(all terms)] = 57.058 (26). (66) 
The total contribution of Set Ill(a) to the muon g — 2, the sum of ([65]) and ([66]) . is thus 

aJ,^°)[Set III(a)(all terms)] = 112.418 (32). (67) 

The total contribution of Set Ill(b) to the muon g — 2, from ([56]) . ([57]) . and ([60]) . is 

a^^°)[Set III(b)(all terms)] = 15.4074 (45). (68) 

The contribution of Set Ill(a) to the muon g — 2 is very large, which is not unexpected. 
In particular, the orders of magnitude of contributions from the dominant (mee) terms of 
Set Ill(ad) and Set Ill(as), as well as the (me) term of Set Ill(b), can be estimated cru dely 
since their leading log{m^/me) term is determined by the renormalization procedure [3], [ssf: 

a(,i°)[Set III(ad)('"^")] ~ 3Kla'f\no loop) - 34, 
4^°) [Set III(a,)(™"")] ~ 3K2,2af\no loop) ~ 34, 
a(,^°)[Set Ill(b) ('"'=)] ~ 3K^a'f\no loop) ~ 7, (69) 
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with 



K2 



2,2 




-ln(m^/me) 



. . . , 



(70) 



and [30 1 



a. 



.f\no loop) = 0.904 979 



(71) 



where no loop means diagrams without closed lepton loops of vacuum-polarization type. 
The factor 3 accounts for the increase in the number of diagrams caused by insertion of 
vacuum-polarization loops. As is expected from (|69|1 and (ITOl) . the values of fl26l) and fHT]) 
are of the same order of magnitude. 
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